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Dynamics of hydrogen atom driven by microwave field of arbitrary po- 
larization is approximated by the discrete mapping. The map describes the 
change of dynamical variables from an aphelion or a perihelion to the next 
one. The results are compared with numerical simulation and previous 
approximations . 

PACS numbers: 05.45. +b, 32.80.Rm 

1. Introduction 

The microwave ionization of Rydberg atoms, a standard topic in the 
discussion of the classical-quantum correspondence in classically chaotic 
systems [1, 2, 3, 4], brings still to us new challenges and surprises. The 
development of experimental techniques allows one now to vary experimen- 
tally the microwave polarization [5, 6]. That in turn enables an efficient 
control of the dynamics and is a challenging problem for theoretical study 
[7, 8, 9, 10]. As it becomes now more and more clear, different quantum 
initial states contribute predominantly to the onset of the ionization when 
the microwave polarization changes [6, 11]. 

The full numerical quantal description of the problem becomes within 
reach of the present days computers. Still simplified models have con- 
tributed a lot to our understanding of the mechanism of excitation of an 
atom and its subsequent ionization especially in term of classical chaos on 
one side and dynamical localization in the quantum approach [2]. Among 
such a simplified description a special place is taken by the Kepler map 
[12, 13, 14, 2] where the full continuous dynamics is approximated by a 
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discrete mapping. The name Kepler map was introduced in [13]. This 
map could be used both to classically estimate the experimental ionization 
threshold as well as to give the quantum mechanical predictions after an 
appropriate quantization [13, 2, 15, 16]. 

The main idea behind the mapping approach is an understanding that 
the electronic motion on the Kepler orbit is the most sensitive to an ex- 
ternal microwave perturbation when the electron passes in the vicinity of 
the nucleus. Such an approach was first applied to a straight line orbit 
(one-dimensional approximation) [12, 13, 14]. This approach was further 
developed and quantized by Casati, Guarneri and Shepelyansky [13, 2]. The 
kicked-like dynamics assumes a Kepler electronic motion between passages 
at perihelions and the kick every time the electron is at perihelion. One 
may then obtain the map describing the dynamics between the successive 
passages at perihelions, or, as proposed in [2], to make the map describing 
the change of dynamical variables from an aphelion to the next one. To 
calculate the kick strength the influence of microwaves over whole Kepler 
period was integrated. This is why the Kepler Map may be seen also as 
the approximation of the mapping generated by Poincare section at radial 
momentum p = [17]. This argument shows that approximation should be 
also valid when the electron orbit stops to be linear. Further approxima- 
tions [2] lead to a map equivalent to a standard Chirikov map [18]. The 
quantum mechanical behavior of the latter is well understood [19] leading 
to an explanation of classical-quantum differences in microwave ionization 
of H atoms in terms of the Anderson (dynamical) localization. 

The authors of [2] extended their analysis to the circular microwave 
polarization. However, soon Nauenberg showed [17] that the dynamical 
variables in the derivation of the Kepler Map are not canonically conjugate. 
He derived the Kepler Map for H atom in circularly polarized microwaves. 

We propose the construction of the Kepler Map for H atom subjected 
in microwaves of arbitrary polarization. This construction allows us to 
answer questions concerning "canonicity" of variables in previous works. We 
analyze the limiting cases of the linear motion of the electron in the same 
linear polarization of the microwave field, and the circular polarization of 
the microwaves. Our new Kepler Map reduces to the previous approaches 
in both cases. We compare the map with numerically integrated equations 
of motion and study validity of the approximation. Further, surprising 
limitations of such an approach are discussed. 



Kepler Map for H atom driven by microwaves. 



3 



2. Construction — one dimensional model 

The Hamiltonian of one dimensional (ID) hydrogen atom in microwave 
fields reads (in atomic units) 

V 2 I 

H = - + Fxcosujt (1) 

2 x 

where (x,p) are the position and momentum variable, F is the amplitude 
of microwaves and u> their frequency. The motion takes place for x > 0. 
Perihelions of the motion described by this Hamiltonian correspond to the 
point x = 0, when the electron bounces from the nucleus. Momentum 
p = dx/dt then is infinite but the whole Hamiltonian H has finite limit for 
every x — ► 0. To eliminate the time dependence of H we go to an extended 
phase space [18] with (H, t) being the second pair of variables. The new 
autonomous Hamiltonian is 

p 2 1 

H = - + Fxcosut- H = 0, (2) 

2 x 

leading to the following equations of motion (£ is the new time of evolution) 

(3) 

S = $ = -!. f = -^ = Fxusmut. 

We evaluate the change of (H,t) variables from a perihelion (x = 0) to 
the next one. These variables are canonically conjugate so the resulting 
mapping should be a canonical transformation. We introduce the eccentric 
anomaly u: cosu = 1 + 2Ex. E is the Kepler energy E = ^ — ^ of the 
electron. The parameter u allows us to describe perihelions (n = 2kir, k 
integer) and aphelions (u = ir + 2/c7r) of the motion of the electron, since 

the momentum is equal to p = ^ — 2Ej^^. We can calculate ^ 

du (-2E) 3 / 2 ( F \ 

' 1 + -=~(1 -cosu) cos ujt) . (4) 



d£, 1- cosu V 2E 2 

If J^- < 1 the derivative ^ may be inverted and the motion of the electron 
can be parametrized by u 

dH Flo (1 — cos u) 2 smut 

*7 = (-2£) 5 /2 1 + ^(1- cosu) cos ' U 

(6) 



2E- 

dt —1 1 — cos u 

du ~ (-2Ef/ 2 1 + ^(1 - cosu) cos ut 
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Integrating these equations to the first order in F/E 2 for < u < 2tt we 
obtain ID Perihelion Kepler Map 



H' 



Flu f 2n 
H + (_2ff/)5/2 / (l-cos U ) 2 sin^rdn, 



(7) 



t' = t- 



2tt 



ujk (-2H>y/ 2 J 



2n r 



2(1 — cos-u) 2 — 3sinu(n — sinn) cosuirdu, (8) 



where r = (u — smu)/ujK + t is the time obtained from equation (6) in- 
tegrated in zero order in F/E 2 with r = t for u = (at perihelion), 
u>k = (-2H) 3 / 2 is the Kepler frequency. The condition of area preservation 
leads to an implicit character of these equations (the right hand side depends 
on H'). This map is generated by the function G (H = t' 



dG 



2tt F f 27T 

G(H' t) = H't + ; + — - / (1 - cos u) 2 cos Urdu 

K ' ' V^2lF {-2H'f/ 2 Jo V ' 



(9) 



The map (7), (8) coincides with Nauenberg ID Kepler Map [17], although 
his map is derived by evaluating the change of Kepler energy E in place of 
H. For the ID dynamics the perihelions always take places at x = 0, where 
E = H. 

The alternative to the perihelion map is the aphelion map. To obtain it 
we must integrate eqs. (5) and (6) for — ir < u < tt with the new function 
r = (u — smu + tt)/ujk + t, since now the condition is r = t for u = —tt. 
The ID Aphelion Kepler Map reads then as follow 



H' = H + 
t' = t- 



Fco 



{-2H'f/ 2 
2vr F 



cosu) 2 smturdu, 



/7T 



CO 



5 cos lot — 3 (u — sin u + tt) sin lot 

UK 



du. 



(10) 



(11) 



This map does not coincide with Nauenberg ID Kepler Map [17] any more. 
Since these equations are related by the condition of area preservation, we 
can concentrate on the first one (10) only After integrating by part and 
some calculation we have 



H' = H + 



-2H' 



LU 



sin ( ujt + tt 

LO K 



4sin7r / sin u sin ( (u— sinn)) du 

JO \UK J 



(12) 
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The remaining integral has nice asymptotic behavior for u ujk, this fact 
serves in [2] to approximate Kepler Map by the standard-like mapping. In 
the high frequencies regime we find 



H' = H + F sin ( ut + vr— ] 

ukJ 

12vr 



-2H' 



U 



2/3 



t' = t F 



2vr 

U K 



(-2H') 7 / 2 



cos 



U 

ut + 2vr 

U K 



+ 



(13) 



(14) 



+ 



7.74vr 



( u ^ 


\ ■ W 1 




sm7r 


V U K J 





As in standard map the dependence of iT' on t is sinusoidal, but the phase 
and the amplitude depend on the energy H'. The change of Kepler energy 
E from a aphelion to the next one is equal to 



w ^ 2.58F / u 
E' = E — - sin ut + vr 



(15) 



but t is not the canonical conjugate variable to E. Disregarding this fact 
one could try to perform a canonical change of variables to have new time 
T = t + tt/ujk, keeping the energy E fixed. Then however, to express the 
mapping in these new variables, first it is necessary to solve an implicit 
equation, which is an analog to eq. (15). 



3. Comparison of one dimensional dynamics 

To verify our approximations we have integrated numerically the equa- 
tions of motion of ID model of hydrogen atom (3) in parabolic variables. 
Along each trajectory we have located all aphelions or perihelions. The 
figure 1 shows the total energies and times (out modulo 2tt) of subsequent 
passages by aphelion generated from the numerical integration, ID Aphelion 
Kepler Map, and the same map in uj S> luk regime. To find the next itera- 
tion of the map we solve numerically implicit equations (10) and (13). For 
H < —0.6 the trajectories, found by means of the map, follow ones coming 
from numerical integration. The approximation breaks down for larger H, 
because in derivation of the map we have approximated the dynamics only 
to first order in . When H approaches to the implicit equation defining 
the map stops to have unique solution, this fact explains the straight hori- 
zontal lines for H > —0.4 on the maps picture. Thus while exact dynamics 
may lead to an unbound diffusion (and subsequently ionization), the map 
cannot be extended to that regime. 
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Fig. 1. Comparison of numerically integrated dynamics of ID hydrogen atom (left 
panel, denoted as a) with the Aphelion Kepler Map (middle panel, denoted as b) 
and the same map in the high frequency limit (right panel, c). The total energy 
H is plotted versus tot modulo 2tt. The amplitude of microwaves F = 0.02 and 
frequency lo = 1.6. 



The similar comparison may be performed for the earlier approximations 
of Kepler dynamics [12, 14, 2]. Figure 2 compares the aphelions found from 
numerically integrated trajectories, the Kepler Map of Casati et al. [2] 
and its high frequency limit version. Here the Kepler energy E is plotted 



a) b) c) 
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Fig. 2. Comparison of numerically integrated dynamics of ID hydrogen atom (left 
panel, denoted as a) with the Kepler Map of Casati et al. [2] (middle panel, 
denoted as b) and the same map in high frequencies limit (right panel, c). The 
Kepler energy E is plotted versus tot modulo 2ir, the Lot variable of the map must be 
shifted by tt for correspondence with time. The amplitude of microwaves F = 0.02 
and frequency lo — 1.6. 

versus tot modulo 2ir. The tot variable for the maps are shifted by 7r to 
restore the correspondence with the true time. The use of E as the map 
variable and the shift of time comes from the definition of the map in [2] . 
No problem with uniqueness of solutions occurs when iterating these maps 
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due to asymptotically (uj S> ujr) explicit character of equations denning 
them. 

A similar picture comparing the ID dynamics with the Kepler Map is 
shown also for higher frequency oj = 3.7. In figure 3 the trajectories gener- 
ated from the maps follow the numerically integrated ones for small H. The 
higher resonances 3:2 and 4:3 {uj : ujk) appears stronger and in a different 
position in the discrete dynamics of the map than in the continuous case. 
Figure 4 shows that for same frequency uj = 3.7 the canonical transforma- 
tion of non canonical variable (E, t) may misplace some characteristic part 
of the phase portrait. The resonance 2:1 is moved by ir in the ujt + n variable 
(the 7r in the variable comes from the definition of the map) . 
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Fig. 3. Same as Fig. 1 but for F = 0.02 and uj = 3.7. 



a) b) c) 
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Fig. 4. Same as Fig. 2 but for F = 0.02 and uj = 3.7. 

We have observed that Perihelion Kepler Map approximates the numer- 
ical dynamics of the model for a somewhat larger range of the energies H. 
The figure 5 shows the Poincare sections generated from perihelions found 
from numerically integrated trajectories and found by means of the Kepler 
Map (as given by the eqs. (7), (8)). The better approximation for larger 
total energy H may come from the fact, that at perihelions the influence of 
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perturbation is maximal, the variables change rapidly. The map, obtained 
by the integration of the influence of the microwaves over the Kepler period, 
should be better if the interval of integration do not cross the perihelion. 
Still the map suffers from the same drawback — it is valid locally only, given 
by implicit equations, which does not have to (and typically do not) yield 
a unique solution far from the initial conditions. Thus no prediction e.g. 
for the ionization threshold may be made basing on the "proper" , canonical 
Kepler map. 
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Fig. 5. Comparison of numerically integrated dynamics of ID hydrogen atom (left 
panel) with the Perihelion Kepler Map (right). The total energy H is plotted versus 
u>t modulo 2tt. The amplitude of microwaves F = 0.02 and frequency to = 3.7. 



4. Construction — two dimensional model 

The two dimensional model of hydrogen atom may be perturbed by the 
microwave field of any polarization. The Hamiltonian of the autonomous 
system in the extended phase space reads 

p 2 I 2 1 

H = 1 x 1- Fr(cos <j) cos cut + a sin <j) sin uj£) — H = 0. (16) 

2 2H r 

We use radial coordinate r, (f> and t as the position variables and p, I, H 
as the corresponding momenta. The parameter a describes polarization of 
the microwaves, it ranges from (linear polarization) to 1 (circular one). 
In analogy to ID model we introduce the eccentric anomaly defined by 
ecosn = 1 + 2Er, where E is the Kepler energy E = ^- + ^ — £ of 

the atom and e is the eccentricity e = a/1 + 2E1 2 . The points of Poincare 
section at p = (perihelions and aphelions) are described by u = kir with 
k integer, since 
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If the perturbation is small F/£ 2 <C 1, we can find ^ and ^ and integrate 
them for < u < 2-k to first order in F/E 



2 



Flo f 2v 

H' = H+- -TTT / (1 — e'cos«) 2 (cos#sinu;T — a sin 9 cos lot) du, (18) 

(-2H') b ' z Jo 

P f2lT 

l' = I + — rrr / (1 — e'costi) 2 (sin0costi;r — a cos 6* sino;r)(iu, (19) 

(—2H') t >/ 2 Jo 

where functions r = (u — e' sinu)/uJK + t, lok = (—2H') 3 ^ 2 and 9 = x( u ) + <ft, 

f u , VT^sinu 

= / l ; jdu , sm X (u) = — , 20 

Jo 1 — e cos v! 1 — e cos u 

the eccentricity e' = VT+ 2H'l' 2 . The equations (18) and (19) are generated 
by the same function 

2tt 

G(H', I', t, (/)) = H't + l'4> + - + 2vr/'+ 

+7 / (1 — e! cos ti) 2 (cos 9 cos lot + a sin sin ojT)du. (21) 

{-2H') b l z Jo 

This generating function defines a 2D Perihelion Kepler Map for hydrogen 
atom perturbed by microwaves of arbitrary polarization 1 . 

Two limiting cases of the transformation generated by (21) may be ver- 
ified. The linear motion of the electron along the axis of linearly polarized 
microwave field corresponds to / = 0, sin</> = and a = 0. One can check 
that the implicit equation (19) for V has the solution V = I = 0, the sin</>' re- 
mains 0, because (ft changes by 2tt, and the generating function (21) reduces 
to the ID Perihelion Kepler Map (9). 

For circularly polarized microwaves a = 1 and the integrated function 
in (21) reduces to (1 — e' cos u) 2 cos(9 — lot). It depends on t and (ft via the 
common variable (ft — cot. We may canonically change the variables to have 
(ft = (ft — cot and t = t, so H = H + loI and 1 = 1. The generating function of 
new variables reads as follows 

- - _ 2-7T 

G(H', I', t, (ft) = H't + I' (ft + ; = = + 2vr/'+ (22) 

-2{H' - loI') 



F f 2v (to ~\ 

= = — / (1 — e' cos u) 2 cos ( y(u) (u — esinu) + 6 I du. 

{-2{H' - Lol')fl 2 Jo V ' \ y ' ujk ' J 

1 The very same procedure is possible in 3D where I 2 should be replace by the total 
angular momentum squared L 2 . 
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H is invariant under the transformation (4). It is the autonomous Hamil- 
tonian of the system in the coordinate frame rotating with the microwave 
frequency. The function (4) generates in variables (I, (ft) Nauenberg Kepler 
Map for circularly polarized microwaves [17]. 



5. Comparison of two dimensional dynamics 

The Poincare sections generated by 2D Perihelion Kepler Map are four 
dimensional. The fast dynamics in (H, t) variables is accompanied by the 
secular motion in (L,cft). Figure 6 shows the (L,<ft) projection of Poincare 
section generated from perihelions found by numerical integration of equa- 
tions of motion for 2D model of hydrogen atom and found from 2D Peri- 
helion Kepler Map. All trajectories start with H = —0.5 and t = 0. Small 
oscillations in trajectories come from the rapid motion in (H, t) variables. 
The shadow of isolated points placed near L = on the Poincare section 




Fig. 6. Comparison of numerically integrated dynamics of 2D hydrogen atom (left 
panel a) with Perihelion Kepler Map (right panel b). The angular momentum L is 
plotted versus the radial angle (ft, this is the projection of Poincare surface of section 
so the trajectories may cross. The amplitude of microwaves F = 0.01, frequency 
uo = 2 and the polarization a — 0.5. 



with numerical trajectories is caused by fact, that the radial angle, at which 
the perihelion occurs, stops to be well defined for very elongated orbits. 
The map does not suffer from this problem. To demonstrate this feature we 
plot (Figure 7) the dependence L{cft) and L(t) in perihelions for a generic 
trajectory crossing L = with same parameters as in Figure 6. The time 
dependence of angular momentum shows that the Kepler Map predicts good 
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values of L at subsequent perihelions. The small error in location of perihe- 
lions in the numerically integrated trajectory results in a much bigger error 
of radial angle, due to a very rapid motion of electron near perihelion, when 
the orbit is almost linear. The time dependence of the angular momentum 
L demonstrates that the Kepler Map evolves a little faster than the true 
dynamics. 




Fig. 7. The dependence of angular momentum L on radial angle <j> (a) and time t 
(b) in perihelions for a generic trajectory. The crosses denote </> values found from 
the numerical integration of equations of motion for 2D model of H atom and the 
solid line is the result of applying 2D Kepler Map. For L w the radial angle at 
perihelion is not well defined, this results in isolated points in the left plot. The 
t dependence does not show a similar problem. The amplitude of microwaves is 
F = 0.01, frequency u) = 2 and the polarization a = 0.5. 



6. Conclusions 

The construction of approximate dynamics of hydrogen atom driven 
by microwave field of arbitrary polarization by means of a canonical map- 
ping has been presented. The approximation - Kepler Map was compared 
with numerically integrated dynamics of perturbed H atom and with pre- 
viously constructed mappings [12, 14, 2]. Although the application of the 
map requires solving an implicit integral equation, which is not faster than 
integration of equations of motion, the map may be used to reduce the 
dimensionality of dynamics and to an approximate quantization. 

On the other hand the maps obtained suffer from a major drawback. 
While being accurate in the vicinity of initial conditions they cease to be so 
when the electron energy increases significantly due to the diffusion. The 
effective microwave amplitude becomes then quite strong and the maps, 
obtained in the first order in F/E 2 , no longer approximate the dynamics. 
Even worse, they do not yield a unique solutions in that regime. This 
behavior may be contrasted with the original Kepler map due to Casati and 
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coworkers [2] which while not canonical and, strictly speaking, incorrect 
yields some estimate for the onset of chaotic diffusion and ionization. Its 
other great advantage is the simplicity of its form. Still, it is our opinion, 
that one should be very cautious while trying to extract any quantitative 
information from this celebrated mapping due to the shown inconsistencies 
in its construction. 
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